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CHAPTER I 


Micromechanics of a Crack-Layer (CL) 

I . Introduction 

In recent years a significant amount of experimental data has 
been accumulated on a process zone surrounding the tip of a propa- 
gating crack [1 ,2, 3, 4, 5, 6]. Process zone is usually defined as the 
area of severely damaged material adjacent to the crack tip. It has 
been shown in works [7,8,9,10] that mophology of damage zone varies 
from one material to the other. In these works ceramics, rocks, 
polymers, and metals were investigated. It was found that damage 
constituting the process zone can reveal itself as microcracking [3, 

II, 12,6] in all of the above materials; martensitic transformation 
in ceramics, metals, and polymers [13,14,15,7], slip lines (i.e., 
shear bands) in metals and polymers [16,17], crazing in polymers, 
etc. 

Despite the difference in morphology of process zones in various 
materials, there are similar features in all of them. For example, 
similar global geometry and similar kinetics of development have 
been observed [18], Theoretical models have been proposed for the 
description of kinetics of a process zone [9,13,16,18,19]. It should 
be noted that traditional fracture mechanics can be considered as 
one of them. Both linear and nonlinear fracture mechanics take into 
consideration the crack tip plastic zone and use the well developed 
techniques of plasticity theory for estimates of its size, shape, 
etc. Damage distribution, however, can be quite different from the 
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FIGURE 3 


one described by models of plastic behavior. The results of work 
[20] clearly demonstrate that some damage patterns do not yield any 
model of plasticity. The theoretical model of a crack layer (CL), 
as opposed to the crack-cut of traditional fracture mechanics, was 
proposed by A. Chudnovsky in 1976 [18], In the CL theory a crack 
with the surrounding damage is considered as a single macroscopic 
entity. The process of crack propagation is described as a 
nucleation, development, and subsequent coalescence of microdefects 
in a crack tip zone. Recently the theory was examined in works 
[20-22] with the goal to observe damage distribution patterns in a 
crack tip zone and to investigate the kinetics of damage develop- 

TV;','- 

ment. This has been done for various materials and different loading 
histories. Figures [1,2,3] taken from [20,21] represent microphoto- 
graphs which illustrate the discussion above. These experimental 
results lead to the schematic representation of a CL on Figure 4. 

ACTIVE ZONE V A 



Figure 4 
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The CL is represented as a crack-cut surrounded by damage of 
density P. Damage density p exceeds the level Po (i.e., p > PQ ) in 
both macrocrack tip active and inert zones; damage density rate p, 
however, is positive in active zone but vanishes in the inert 
zone Vp 

Instead of a detailed description of damage density p and its 
evolution, the CL theory operates with integral characteristics of 
the damage zone. These integral characteristics are: the length, 
the width and the shape of the active zone. The main crack length 
and the curvature of the crack trajectory also enter as geometric 
characteristics of a CL (we talk of crack trajectory rather than 
crack surface when considering the crack within the framework of a 
plane problem). On the basis of general principles of irreversible 
thermodynamics, the authors of works [18,23,24] have introduced CL 
driving forces which are reciprocal to geometric parameters of a CL. 
the detailed description of a CL geometry and kinetics can be found 
in [25]. This work introduces four CL driving forces which appear 
to be reciprocals to 1. CL length, 2. the curvature of CL trajec- 
tory, 3. active zone area, and 4. active zone shape (active zone 
length to width ratio). The first three CL driving forces are 
represented as linear combinations of the well known path-indepen- 
dent integrals J, L and M of fracture mechanics and, the integral N, 
[25] which is not path-independent in general. 
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The Integrals J, L, and M appeared in the formulation of conser- 
vation laws of el astostatics , (i.e. J = 0. L = 0, and M = 0 on any 
closed contour without singularities inside of it, (263). These 
conservation laws have been derived by means of Noether’s theorem 
from the principle of minimum of strain energy. The conservation 
laws J * 0, L = 0, and M = 0 result from invariance of strain energy 
functional with respect to the group of displacements, rotations, 
and infinitesimal isotropic expansions, respectively. Naturally, 
the first one holds for homogeneous material only, the second holds 
only for isotropic and the third one for linear medium. The conser- 
vation law involving J-integral appeared for the first time in work 
[273 of 1951 by Eshelby. In terms of the J-integral Eshelby express- 
ed the force acting on a singularity in an elastic body. Later, the 
J-integral was rederived by Sanders [283 in 1960, Cherepanov in 1967 
[293, and Rice In 1968 [303 in connection with a problem of energy 
release rate in a quasi static crack propagation process. The tan- 
gential to the crack component of vector J was found to be precisely 
the energy release rate. After the discovery of the other two 
conservation laws, L and M were interpreted as energy release r^tes 
with respect to cavity rotation and cavity expansion, respectively 
[313. In works [23,24,253, by A. Chudnovsky et al , the path-inde- 
pendent integrals appeared as parts of thermodynamic crack driving 
forces. 

For evaluation of J, 1* M, N integrals within the framework of 
Cl theory the field around the CL must be known. In order to evalu- 
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ate the stress field a mathematical model of the microstructure of 
the CL must be introduced. Such a model has been proposed in [32]. 
This model considers the micromechanics of the CL, i.e. the 
interaction of a macrocrack (or main crack) with an array of 
microdefects in a close vicinity of a macrocrack tip. The work also 
outlined the method for evaluation of a stress field around the CL. 
The present work develops the method of [32]. 

The model of work [32] treats a CL as a crack-cut in linear 
elastic medium surrounded by microcracks in the crack tip zone. The 
microcrack array surrounding the macrocrack tip consists of a field 
of randomly distributed small cracks (small in comparison to the 
macrocrack) of random lengths and orientations. The problem is to 
find the elastic stress field resulting from interaction of an array 
with the main crack subjected to external tractions. The main 
crack, with the microcrack array surrounding it, has been considered 
within the framework of the plane problem of elasticity. Also, the 
assumptions of small scale microcracking were supposed to hold [32]. 
This means that the microcrack array occupies a small area in com- 
parison with the main crack tip. Under this assumption the macro- 
crack asymptotic stress field appears to be dominating and it is 
defined by the stress intensity factor K° only. - , 

For the sake of simplicity, the mode I loading conditions have 
been assumed to hold. This explains the notation for intensity 
factor K°, superscript "0" refers to the main crack. 

The method of potentials, has been chosen as a means of solution 
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of the problem. The reasons for that selection are: the possibility 
of generalization of the method for 30 problems, and the convenience 
of the method for statistical purposes. 

The method of potentials gives the elastic stress field in a 
form of integrals of the potential density multiplied by Green's 
function and, therefore, the solution explicitly depends on the 
microcrack array configuration. This form of a stress field solu- 
tion permits a relatively easy statistical averaging procedure. 

2 . Mathematical Formulation of the Problem 

The macrocrack (or main crack) interacting with an array of 
microcracks under the assumptions of small scale microcracking 
(definition follows) is being considered. The two-dimensional 
linear elastic solid contains a macrocrack of length "2£ o " and an 
adjacent array of N rectilinear microcracks of lengths "2 % each, 
with n. as a unit normal vector to the i-th microcrack, and as a 
position vector on the i-th microcrack, i = 1, 2. .... N. The 
elastic plane is under mode I tensile loading with respect to the 
main crack. 

The stress field a(x) can be represented as a superposition 

N 

0(x) = o^ + .qCx) + E =1 S^*) l 1 * 1 ) 

where o°° is the stress field due to remotely applied loads in the 
absence of cracks, a and 0 .(x) are the stress fields generated by 
the main crack and by the i-th microcrack, respectively. More 
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exactly, is the stress field in an infinite solid containing one 

crack (references will be made to the microcrack which is just i- 

th microcrack) with faces loaded by tractions n. [o°° + G(x.) + 

N : ~ 1 ' 

|^?i ) ] where o(X^) and a^(X^) are actual stresses generated 

by the main crack and the k-th microcrack along the line of . 

In the vicinity of the microcrack tip, stresses can be neg- 
lected compared to the tip-dominated field 


2(x) 


K eff 4>[6(x) ] 
^ /2nr (x) 


rff 

Where denotes the stress intensity factor for the macrocrack 
tip with the effect of microcracks taken into account, r and 9 
denote the position - vector and polar angle in the main crack tip 
corrdinate system. > The small scale model is defined by the condi- 
tion a m << d(x). Thus, the asymtotic stress field in the vicinity 

of the main : crack tip can be represented in the form. 

d>( 6 (x) ] ff N 

o(x) - -=== — Kj_ + - i £ 1 2 i <?) (1.2) 

/2nr (x) 1 1 


The technique of double layer potentials, with a potential 
density as a crack opening displacement (i.e. displacement discon- 
tinuity) will be used [33]. The displacement field may be represent- 
ed by means of double layer potential density as follows: 



where b(^) is the crack opening displacement, and x) is the 
second Green's tensor of elasticity for plane stress conditions. 
(Here and below through the whole work d£ must be understood as the 
increment along the crack.) The tensor 4- ( c, ,x 1 constructed in 
Appendix I may be written as follows: 


J(C,x) 


1 + v 

= ,~7T [d - 2v)(n R - Rn - n RE) 

4 n R x -X — - 


n x JR 

- 2 — y-'RR] 


(1.4) 


where R = 5 - x, v - is Poisson's ratio, E - is the second rank unit 
tensor, and the factor 1 + v is to be substituted by the factor of 
1/(1 - v) for plane strain conditions. 

Differentiating the displacement field (1.3) and taking the 
symmetrical part of a tensor gradient, the strain tensor can be 
obtained. Application of Hook's law to the strain tensor results in 
the stress tensor o(x): 


o(x) 



b(p-4>(§,x )d £ 


Z 


0.5) 


Where J x is the stress operator transforming the displacement 
field u(x) into a stress field (subscript x indicates that differen- 
tiation in T x is performed with rp spect to x). Thus, the stress 
field ( 1 . 2 ) can be represented in the form 


o(x) 


$[9(x)l K eff + 
/27Tr(x) 1 


N 

ill 


T 

-x 


j b.(C)-$(C, X)d§ 


( 1 . 6 ) 
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where N unknown functions b.(Z) are to be determined from N vectori- 
al integral equations expressing boundary conditions on the micro- 
cracks The equilibrium equations are automatically satisfied 
for both terms in the superposition formula (1.6). The first term 
yields equilibrium equations because of the properties of the 
asymtotic crack solution, the second - because of the properties of 
the second Green’s tensor $(£, x). 

The faces of the microcracks must be traction free: 


2 i ( S < Si ) + K ii v 


N 

A 

K^i 


K 


b K (§)*H^.Xi> d ^ + T x b i (0*$(§»^ i ) d ^ " 0 

£ i 0.7), 


for all x e for each i. 

The last integral in the braces is to be understood in the 
principal value sense. It should be noted that the last integral in 
(1.7) becomes divergent if stress operator is moved under the inte- 
gral sign and applied directly to the Green's tensor <t>( 5 , x). It 
can be shown, however, that the limiting value of the integral is 
given by the following regularization: 


8,1111 T 
~ > 

x-*x . 

~ - 1 


J b(§)'*U. 

i 


x)d£ 



[b(£) - b(x)]‘T x $(§, x)d£ 


( 1 . 8 ) 


where the integral on the right must be understood in the principle 
value sense [34]. 

ef f 

Expression (1 . 8 ) contains one more unknown - .An additional 
equation reflecting the impact of the microcrack array on the main 
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(1 .9) 


crack may be written in a form 

i r n 1 

1 ro £ + x N 

= K° + ~?= I— n(x) I o.(x)*n(x)dx 

o 1 l -x" i-1 ’ 

-l r ° 

o 

in the small scale microcracking model. 

The last equation is an exact substitution for the boundary con- 
dition on the main crack. 

Thus, the system of 2N + 1 scalar equations (1.7,1 .9) represents 
the formulation of the main crack-microcrack array interaction 
problem. In the following the system of equations (1,7 ,1.9) is 
solved for two and three crack interaction problems. 

In Chapter II, three particular problems are solved under the 
assumption of macrocrack dominating stress field to be piecewise 
constant on each microcrack. This is an approximate solution of the 
basic system of equations which can be justified for ratios 'f^/t << 

1 i = 1 , 2, ... N) . 

In Chapter III, the higher order approximations ahe considered 
for two crack interaction problems. It is shown that the method can 
be extended to higher approximations and become exact in the limit- 
ing case.' 

In the final Chapter IV, the expression for stress field in a 
general problem of interaction of a macrocrack with the microcrack 
array of arbitrary configuration is obtained. It is shown that the 
resulting stress field can be fully characterized by asmytotic 

stress field of the macrocrack c (x), (in absence of the microcrack 

- n ~ 
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array) the values of its derivatives in the directions of micro- 
cracks evaluated at the centers of microcracks, and the second 
Green's tensor $(5, x) constructed in the appendix I. 

The obtained elastic stress field solution has been used for 
evaluation of the J-integral for the CL. 


14 



CHAPTFR I! 


Piecewise Constant. Approxima t ion i n Two 
and Three Crack Interac t ion Problems 

1 . Two Track Interact i on Problem ( t.wn c ol inner crack s) 

In + his section. t. h e problem of plastic interaction of a macro- 
crack with one microcrack located on the samp lino is considered. 
Thp problem is solved based on t.he formulation of the previous chap- 
ter (i.e. planp stress, small scale model, mode I loading condi- 
tions) when the dominating stress field o(x) on t.he microcrack a. is 
approximated by a constant. 



Figure 5 
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Under modp I external loading, t.hp dominating st.rpss field of thp 
main crack of lpngth ?z^ is given by 


§00 


K eff $ [ e Cx) ] 
^ /27Tr (x) 


( 2 . 1 ) 


Thp stress field duo to t.hp microcrack Z is givpn hy (1.6): 

c+£ 

0 2( x) = T x J b(0'*a. x)dC (1 . 6 ) 

c-2. 

where b(£) is thp microcrack opening di spl acempnt. . (i.p., donhle 
layer potpntial donsity on z ) . In thp case of a constant, approxima- 
tion. microcrack Z is found to hp pmhpdded into a uniform stress 
field. The COO h(£) in a uniform stress fipld is known to he ellip- 
tic r 35 1 and may be written for mode I conditions: 

AO 

= — n(c) *o(c) e (Q (2.2) 


where F-is Young's modulus, n(c)-is unit normal vpctor to thp crack 

£ at its cpnter, o(c)-is the resulting stress tensor at the micro- 

2 

crack center, and e(0 =. 1 - is the elliptic crack opening 

JT 

(n(c).o(c) is a traction vector at point c). 

The boundary conditions (1.7), in view of (2.2), gives 


° 22 (x) = 


° 22 (x) = 


.eff 


/2tt()c + 6) 


x£(c-£, c+£) 


(2.3) 
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where superscript ‘V refers t.n the microcrack. The equation (2.3) 

takps scalar form because of thp symmetry of the svst.pm of cracks. 

Thp strpss componpnt appears to bo proportional t.o K™ on 

thp microcrack (noto that 0 Mx) = °oo(c) whpn x (c - £. c + £) 1 . 

22 ?? 

Therefore, h{£), given by (2.2), also hpcnmps proportional t.o 
off 

. Thus, thp last, pquation (1.9) for determination of p f foctivp 

stress intensity factor becomes linear algebraic pquation with 
off 

respect, to . Fquation (1.9) may be written as follows: 


- K? + 

1 i j%r 




a o + x £ 

0-_ (x)dx 

£„ - x 22 


(2.4) 


In the last equation, o ^(x^ must hp calculated from (1.6) with 
(2.2) as a double layer potpntial density. 

Tho displacement vector u(x) is given hv 


, , At 
y(>0 = — 

E 


] e(£)n'(£) *0(0 x)d^ 

£ 


( 2 . 6 ) 


where thp product of tensors in the integrand may he writt.pn in 
index notation as n a U) *° a gU^ U,x) with a, 0 . y - 1,2 (summa- 
tion by the repeated subscripts is implied). In the coordinate 

system of Figure 5 the integrand becomes: p(f; ) . ( a ^ + <5> 2a^‘ 

£ 

For evaluation of the stress component (x) only the derivatives 
3u-i 3Uo 

r— - and t-* need to he calculated due to the Hook's law: 


22 


. 9u 9u_ 

. V , 1 , 2 S , 2, 

Pl l - 2v^9X 1 9X 2 ; 9X 2 J 


( 2 . 6 ) 
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au l 

dX~ 




A£ 

E 


A£ 

E 


c + £ 
’ 

. 

c- a 
c + £ 


c- £ 


e(0 (a 


3$ 

21 W, 


11 


+ o 


3<J> 

22 W, 


21 


)d£ 


3<t 


e(Q (o 


12 


3<5> 


+ o. 


22 


21 3X„ 22 3X 


)d£ 


(2.7) 


where 


* n (Z>'V 




1 + v 
2 

A m R 


[(1 - 2v)X + 2 -|(q " V 1 


R 


, - 2v> < 5l ' Xj) - 2 - X 1 >1 


A n R 


4> 12 (C,X) 


1+ v 

2 

AtiR 


((1 - 2v) <£ - X ) - 2 (Cj ~ VI 


( 2 . 8 ) 


1 + v 

$> ~(£>X) * - 2 

22 ' A tt R 


[(1 - 2v) X„ + 2 —~] 
Z R 


7 . 7 7 

in the chosen coordinate system, with R * (C ^ - x^) + (the 

subscripts "V" and "2" rpfer to the x and y axis, respectively). 
Substitution of (2.8) and (7.7) into (2.6), results in 


° 22 (X) = 


f \ c + £ 
o 22 ( c ) r 


n 


c - £ 


/-^ + 2c£ - (c ^ - £ 1 2 ) 


(C - X)‘ 


d£ (2*9) 


where x and 5 are coordinates on the horizontal axis. 

Integration gives: (The integral in (2.9) is evaluated in Appen- 
dix II) 
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0 22 (X) ■= 0 22 (C) 0- ■ - 1) 

V X-c ; 


( 2 . 10 ) 


wherp x - c >_i. 

Substituting (2.10) int.n (2.4) we obtain as follows: 


o„(C) 

K? ff - K? + 4=- 

1 t /_ n 


(C) f “/* 0+ X ! 

=" i l o - x 

0 -r x " c 

0 


- l)dX 


(2.4a) 


Tbp second t.prm in (2.4a) roprosonts an increment. AK of stress 
intensity factor due to the presence of the microcrack. 


°22 (C) f 0 * 0 + * 
Mr > l o - x 


- l)dX 


(2.11) 


The term decreases with the increase of c (when x - c >> i.) and 


tends to infinity when c -*■ 2, + l. i.e., when 6 0 (the distance 

0 


between the adjacent crack tip tends t.o zero). 

off 1 

The combination of (2.3) (o 0 -(c) = KT 

?? 1 vOTFTT 


) and (2.4a) 


P T T 

results in linear alqpbraic equation for K-| : 


Kf" - K° + 


i . ’ K i + v <i (6 ^> 


( 2 . 12 ) 


where 


q(6/£) 


+ 6 ') 


r l + X 

1 - X 


f'l -(=M‘ 

X-C 


(2.13) 


and £* = 2 ,/fc . 6* = 6/£ o , = 0 arr> nondimensi onal i zed coordi- 

nates with respect, to the main crack length i From dimensional 
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considerations it follows that q = q(6/«.) depends on ratio of the 
only two characteristic lengths of the problem and i) . 

The solution of (2.1?) is obvious 


,eff 


K 1 

1 - q 


(7.14) 


Th 


/eff .„0 


e graph of /K-j is given in Figures 6 & 7. The graph in 


eff lu o 


Figure 7 is presented in order to illustrate the behavior of /K^ 


for smal 1 6 /£ . 

The oht.ained result indicates that effective stress intensity 

eff 0 

factor K-j increases from to infinity when the distance between 
two cracks tends to zero. 

General superposition formula for stresses (1 .?) can be written 
now as follows: o(X) = 6(X) + 0 £ (X) = 



/ <He(x)i 
i /2-nr (X) 


+ 


1 

/2ti(J, + 6) 


c + fc 

42, f 

T " (c) Tx j e( °t 

c - t 


(?.15) 


Formula (2.15) with given by (?.14) represents the approximate 

solution of our problem { pi ecewi se-constant approximation of the 

resulting stress field o( x) . ; 

■ The first term in the solution (2.15) (K?^ ) represents 

/2iTr(x) 

the dominating stress field of the main crack. The second term 
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Figure 6 




o.oi ai 
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Figure 7 
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Figurp R 

represents thp stress field of the microcrack with the magnitude of 
microcrack opening displacement determined hy the magnitude 
of dominating stress field at thp center of the microcrack .Figure R 
gives qualitative graph of the stress distribution following from 
formul a (2 .15) . 

One more comment should he made with respect, t.n the solution of 
the system of equations (2.3), {2 .4) . This system has been reduced 
to single equation (2.12) for .K^ . (Of course, similar equations 
can be written for o^(c)). In. the problem under consideration, 
i.e.. problem of interaction of two cracks, t.hp obvious solution of 
(2.12) is given by (2.14). However, in more complicated situations, 
i.e. in many cracks interaction problems, it may h P useful to try an 
approximate methods for solution of a system of equations corre- 
sponding to (2.3), (2.4). It is easy to see the meaning of certain 
approximations in the case of simple problem under consideration. 

The [321 work suggested an iterative procedure as an alternative 
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to an exact solution. In our problem, equation (2.12) can be solved 

by means of an iterative process. If one takes as a zero approxi- 
eff 

mation for in (2.12), then an iterative process gives geometric 

series 

K I ff = K i (1 + * + q 2 + q 3 + ...) (2.16) 

and the sum of this series for |q|< 1 coincides with (2.14). On the 
other hand, formula (2.14) is meaningful only for |q|< 1, therefore 
(2.14) and (2.16) are equally valid. 

The sequence of iterated terms has clear physical meaning: the 
first term gives the intensity factor of a main crack k° in the 
absence of a microcrack, the' second term accounts for first order 
interaction, i.e.. microcrack, being imbedded in the field of main 
crack K°, gives the correction to of a magnitude k°q. The .third 
term accounts for double interactions, etc. 

Substituting (2.16) into (2.15) we obtain the formula each term 
of which can be interpreted by means of diagram in Figure 9. The 
first term of the sum represents the stress field of the main crack 
under the loading in the absence of the microcrack, the second 
term represents the stress field of the microcrack embedded into the 
main crack field, the third term represents the correction’ to the 
main crack stress field resulting from presence of the microcrack 
which is embedded into the stress field of the main crack, the 
fourth term represents the correction to the microcrack stress field 
resulting from the correction to the main crack field, etc. The 
solution 
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different senses: the approximation of a stress field o(x) by a 
constant on the microcrack, and an approximation due to a number of 
physical interactions between the cracks taken into account. 

2 . Two Crack Interaction Problem (two parallel cracks) 

In this section the problem of elastic interaction of two paral- 
lel cracks is considered, (see Figure 10) i.e., macrocrack of a 
length "2)i 0 " and a microcrack of length "2a". 
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Figure 10 

Constant approximation of the resulting stress field o(x) on the 
microcrack has been used again. For simplicity purposes instead of 
a double layer potentials technique we use the asymptotic stress 
field solution for the microcrack. 

The purpose of this section is to show that the method works for 
the problem which has not been solved so far (to our knowledge). 
The solution has been obtained under the simplest possible assump- 
tions as an illustration of the method. 

The system of cracks in Figure 10 is assumed to be under mode I 
tensile loading. The small scale microcracking assumption also 
holds. The dominating stress field is given by (2.1). The 022 stress 

component of the microcrack field in. r (x), e (x) coordinate system 

£ £ 

may be written as follows: 
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£ , s *22 [ M?>] 

o 22 00 - Kj - g, i - 

/2Tir £ (X) 


+ K, 


v^ur (X) 


(2.17) 


where the first term represents the contribution of the mode I 
asymptotic stress field, and the second term accounts for the mode 
II, (The resulting stress field a(x) in this problem contains both 
modes of loading, of course.) The stress intensity factors K, and 
K,, are given by 


K 1 °22^ C,6 2^ >/ ^ * K 2 “ °2i (C»<$2^ 
and 


(2.18) 


6, 

^22^£^ = d + s * n ~ 


36 f 


2 Sin T K ^21 (e H> 


V - e £ • - 30 { 

COS — sin — sin • — 

2 2 2 


The resulting stress field a(x) acts as an external field ap- 
plied to the microcrack and under the constant stress approximation, 
equations (2.181 hold. 

The boundary conditions on the microcrack (1.7) take form 


° 22 (C - e 2 > * K^ ff 


'2 nr 


°2lC c .« 2 ) - K* ff W 


/27ir 


(2.19) 


” 1 ^ 2 2 

where 6 o = tan and r 0 = ($•,_ + £) 2 + 6^ are coordinates of 

the center of the microcrack in the main crack tip coordinate 
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system. The equation (1.9) for effective stress intensity factor 
appears as follows: 


K*" - K° + o 22 (C,6 2 > 


\/I f 2 + V c -V|/r f i 


( 2 . 20 ) 


where 


fi. . 


-fi 


l + X 
o 

fi - X 
o 


4» 21 l e ? (x) l 


/2-nr £ (X) 


dx 


( 2 . 21 ) 



-fi, 


22 y dx 
>'2r,r jl (X) 


where relations e f . = 0. (x) and r = rj(x) are given by (see Figure 

ji K, Jm ■ 


1°] . r t (X) - (t 0 + 6 X -X) 2 4- « 2 2 , and 0 £ (X) - sin j- 


"l 

7x7 


' v ♦oofe.) 4> ?1 (9 0 ) 

The factors ^(C,^) ~ — and ° 91 (C,<S 9 )' = £ 


/2xr ^ ^ /2iTr^ 


in (2.20) 


characterize the o^x) component of stress on the main crack.' 

Weight functions f^ and f^ depend on the microcrack location and 

rff 

orientation. The dependence of K-| upon the length of the micro- 
crack is given by the -factor of /-L . Equations (2.19), (2.20) 
constitute a system of three linear algebraic equations with three 

off 

unknowns - c^C, 6 2^’ a 21^ r '’ '**?) and K 1 * Its sol , ution y ie1ds 
formula (2.14), where 
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q(6 1 /£, 6 2 /H) 




lf 2 *22 (6 o> +f l *21 (0 o> 


(2.22) 


and stress components are given by (2.19). 

General superposition formula (1.2) with the help of (2.18) and 
(2.19), gives the asymptotic stress field in the vicinity of the 
microcrack in the form: 


o(X) 



4>[6(X)] 

{ + 

v^Trr (x) 


1 2r *t <l> 22 (0 o ) 


<-t e iL (x ) ] 

/2mr fc (X) 


+ 4> 21 C0 O ) 


! [9 £ (x)] 


r2irr^(X) 




The presence of mode II loading in the resulting stress field 
ef f 

o(x) gives rise to K-| ^ in this problem. The shear mode contribution 

is represented by the last terms in formulas (2.22) and (2.23). 

ef f 

In order to make explicit the dependence of K-j on the para- 
meters of a crack system (Figure 10) the function q(6^/£, 6^/ £) must 
be evaluated. This function, in fact, depends on ratios of S^/n, 
6^/1 only, as follows from dimensional analysis. 

Let us consider the expression. 


a 

■ 0 


'2 m r 


__1 

/2mr 


o -£ 


ISL + X 
o 

£ - X 

o 


/2mr ^(X) 


dx 


(2.24) 


The integrand in (2.24) has a singularity at x'= £ q (at the right 
main crack tip) of a square n 


restricted by the inequality 0 _< 


^22[%00] 

^nr t (X)" 



TT 


i 

X 

■ 


1 .3 
/2mr, 


which may be 


narrowed for | e ^ | < m (this is our case because the interaction of 


only left microcrack tip with the main crack has been taken into 


2.23) 
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account). Thus, by the mean value theorem (2.24) may be represent- 
ed as follows: 



(2.25) 


(note that 




-it Z o 


- dx = £ tt) where 0 < a < 1.3. 
X o - - 


reasoning i$ applicable to the second term q(6^/x 


The same line of 
5 2 /0 i n formula 


( 2 . 22 ), 


h m i_ 

/2-nr /2nr 



* 21 ie 8 <X)] BS o 

dx = 

/2nr (X) 2r 

x. o 


(2.26) 


where -0.4 < £ < 0.2. (interval for 6 can be narrowed also). 

Using (2.25), (2.26), and (2.22) q = q( 6 -j /«. , 6-,/fc) takes the form 


In £ ' 

q ( V £ ’ 6 2 H) 1 ~T~ * [ai 22 (0 o ) + 6< ^21 ( 6 o ) 1 


(2.27) 


Introducing r Q = r Q (£, 6^, 6^) into the last expression and 
dividing both numerator and denominator by z Q results in 


q (6_ /£, 6 _/£) = i [out, (9 ) + 64*,. (6 )] 

i t i it o 21o 


V 


(6 X + vy + 6 2 


~ 2 (2.27a) 


■ i 

where z' = a /£ Q , s-j = 5i / £ 0 , and = S 2 /S. 0 are nondimensional ized 
parameters. 
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Direct calculations show that the factor of ( © ) + 
0 *21^ 0 o^ is P° sitive for any location of the microcrack, which means 
that effective stress intensity factor K® ff always increases because 
of the presence of the microcrack. 

In the case of both, 6 ^ 0 formula ( 2 . 22 ) for q may be 
shown to reduce to the corresponding formula of the previous sec- 
tion. Thus, the effective stress intensity factor tends to infinity 
when the distance between micro and macrocrack vanishes. 

Qualitatively the behavior of effective stress intensity factor 
in this problem is described by (2.14) with q given by (2.27) or 
(2.27a). 

3. Three Crack Interaction Problem 


r 



Figure 11 
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The purpose of thissection is to analyze the problem of inter- 
action of a macrocrack with two symmetrically located microcracks. 
Figure 11, under the assumptions of the previous two sections. 
General i zation in both directions nonparallel system of cracks and 
nonsymmetrical location of microcracks is possible but calculations 
become considerably more complicated. 

Two distinct cases may occur, one when the microcrack array 
amplifies the dominating field (and this has happened in two pre- 
vious problems), the other when the microcrack array diminishes the 
main crack field. The latter is called the "shielding" effect. The 
simplest situation for "shielding" to occur is when two parallel 

microcracks are placed at a distance c = sl o (see Figure 11). The 

eff 

goal of this section is to estimate the effect of K-j reduction 
when microcracks have been placed at the indicated position, and to 
show that in case of c >> £ Q and h reasonably small two parallel 
microcracks act in a way similar to one microcrack in section K 

There is no shear mode in the resulting stress field because of 
the symmetry of the problem. That is why the system of equation's 
(1.7) (boundary conditions on microcracks) the form of two scalar 
equations 


° 22 <C - 5 > + F 21 °22 < C -- 5 > * 


- K 


eff 2^ 

^2vr (C , b) • 


(2.28) 


° 22 (C> - 5 > + F 12 ° 22 (C> * 


- K 


eff 


<f> 22 [6(c > ~|>] 
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The first term in the first equation (2.28) represents the 
normal stress component at the center of microcrack the second 
term represents the stress component ^ which microcrack exerts 
on £ i , with as a scalar influence function, and finally, the 
right hand term represents the dominating stress field at the center 
of microcrack The equation is formulated at the center of micro- 
crack because of piecewise constant approximation assumption. 
The second equation has been formulated for the microcrack and 
because of symmetry it is identical with the first one ( o^(c,^) = 
°22 ^ C ’ " | )• The influence function F^ = -j * F (because of 
symmetry) appears from (1.6) with (2.2) as a double layer potential 
density (because the unknown field o(x) is assumed to be constant on 
a microcrack) . It is defined by the expression: 


F(£, x)o = 


4_ 

E 


l 


n(X) n(X)T x j 

-k 


e(C) n(O-o(O-0(c. x)dT 


(2.29) 


where ^ ^ = £ * and n(*) is a unit normal vector at point x in 
the direction of x^ axis. 

The expession above represents the influence function of a crack 
at any point x. The function F( , x*) in a coordinate system of 
the center of the microcack takes the form 


l 

F(C, x’)'= i J e(0 

-9 



- (C - Xj) 4 - (C - x’) 2 xj 2 
HC ~ x’) 2 + x^ 2 ] 3 




(2.30) 


where position vector x' has the components x^ and x^. (The integra- 


32 



tion procedure is described in Appendix II). 

In addition to the system of equations (2.28), there is an equa- 

A f f 

tion for the effective stress intensity factor (1.9). For this 

problem it can be written in the form 



+ 


2 




(X)dx 


(2.31) 


o 

where the factor of 2 in the second term appears due to the presence 
of two symmetrically located microcracks, and o^^(x) with help of 
(1.6) and (2.2) may be represented as 


°22 (X) = T n(X) n(X) a(C * - (C> -x J (2.32) 


where o(c, ^ ) = a(c, - !j) , and n(c, - jjj) = n(c, - ^) , and because of 

that (2.32) holds for both microcracks. Thus, three equations 

(2.28), (2.31), in view of (2.32), represent a system of linear 

h 

algebraic equations for determination of three unknowns °22 (c *2 ) ' 

I -TC 

°22^ C ’ and * (Because of symmetry, in fact, there are only 
two equations for two unknowns (^(c, and 


'? ff 


^ 22 [0(C,2 )] 

Substitution of o..(c, s) = - - t 1 ’ , - iv, 

22 2 V + F(z,0,h) /2nr(C,0) 


from (2.28) 


into (2.32) with the subsequent . substitution of the latter into 

off 

(2.31), results in the equation for K-| with the solution in the 
usual form (2.14) with, . 
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4>„[e(c, g)] { 


q = - 


•nH [1 + F ( 8,0,h) ] 


/2Trr(C,S) 


8. + X 

o 

8 - X 

o 


F(8, X)dx 


-8 


(2.33) 


Influence function F(8, x) in (2.33) is defined by (2.29) and must 
be evaluated in the coordinate system of the main crack (Figure 11). 
This function may be obtained from (2.30) also by coordinate trans- 
formation x^ = x^ - c, = *2 “ £ and integrating from (c - 8 ) to 
(c + 8) . 

In order to analyze the behavior of the effective stress inten- 
sity factor K-j the influence function F(8, x) must be evaluated. 
The integral in formula (2.30) can be evaluated exactly but it is 
enough for our puposes to estimate it with help of the mean value 
theorem: (Formula (2.34) is obtained from (2.30) in Appendix II). 


f ( 8, X') = 





(X -J)/X, 

t 2 '2 

1 + (X x - 8 ) 


(x' + 8) /X ' i 

1 — ^ ) + 

i 2 ' 2 

1 + (Xj + 8) /X 2 


' 3 ' 3 

(X L - 8) J /X 2 J 

[1 + (xj - 8) 2 /X 2 2 ] 2 


' 3 ' 3 

(X L +.8 ) J /X 2 J . 

[i + <x a + 8 ) /x 2 z r 


(2.34) 

)} 


where e Q = /\ _ (^L-L-iT) is the elliptic crack opening at some point 
x Q e(c - 8, c + 8). The influence function F(8, x' ) given by (2.34) 
can be considerably simplified for both large and small 


F(8, x ) = 6 ± 


(2.35) 
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(2.36) 


for large X£, i.e. (X^ -a ) /X^ << 1 and (X^ +Jt)/X 2 << 1. 

F<£, X ) z - £(— r— r~ — ) 

X 1 - l Xj + l 

for small X^, i.e. (X-j -O/X^ >> 1 and (X^ + >> 1. In formulas 

(2.35) and (2.36) it was assumed that e = 1. Substituting formulas 
(2.35) and (2.36) into (2.33) the estimates of K®^/K° for small and 
large £/h can be obtained. It follows from (2.35) that 


F(i, o, h) r 6 ^ 


(2.37) 


for the microcracks being wide apart from the macrocrack tip, i.e. 
£<<h. 

Substituting (2.37) into (2.33), and taking into account that 
the second of (2.37) holds only for - h < X-j j< ft , we obtain 


^ Tt 


6. 

2 

1 + 6(*/h) 


, i ^ 
W 


for small l/b 


(2.38) 


The last expession gives the uppe estimate of q=q(£/h) because of 

the conditions e=l , and % - h < X. < £ . 

o I — o 

Analogously, for the microcracks being close to the macrocrack 
tip, i.e. £ >> h, the formula (2.36) gives 

F(£, 0, h)-*2 for small h (2.39). 

(The last expession represents the limiting value of the influence 
function (2.29) which is 2e Q for the approximation (2.34)). 
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Substituting (2.36) and (2.39) into (2.33) and evaluating the 
integral (see Appendix II) we obtain 


q 


r 





2 + i 

■ — - i _■ i • 

2 


Z 

h 


for large Z/h 


(2.40) 


The important feature of (2.38) and (2.40) is that both of them 
give negative value for q which, in turn, gives the reduction of the 



h/2 

Figure 12 • • ; v 

effective stress intensity factor K^. The graph of vs. ^ 

(i.e., the microcrack length over the vertical distance between the 

macro and microcrack, see figure (11)) is represented in figure 02). 

The upper and the lower portions of the curve were calculated 

using (2.38) and (2.40), respectively. The effective stress inten- 

o o 

sity factor is always less than K-j. It varies from K-j for 
small 4 /h to zero for large 4 /h. 

The effect of reduction of the effective stress intensity factor 

pf f 

K-l at small h is similar to the one obtained on Dagdal e-Barenbl att 
model. The effect of vanishing of the stress intensity factor on 
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Dagdal e-Barenbl att model has been obtained by means of introducing 
interaction forces between the crack faces in a crack tip zone. In 
our model the effect appears as a consequence of microcracking in a 
crack tip zone. 

In the case when the system of microcracks is being moved far 

from the main crack tip in a horizontal direction, i.e., c >> i , 

o 

formula (2.30) gives negative value for F(« , x) , and q in (2.33) 

eft 

becomes positive. Thus, the effective stress intensity factor K-j 
increases in this case, just as in the case of one microcrack on 

the same line with the main crack. More similarity can be noticed 

by increasing the distance between the macro and microcrack centers 

at the given h. In this case the influence function FU » x.) in 

(2,33) may be reduced to the corresponding expression in the problem 

of Section I. , In terms of function F(? , x'), formula (2.30) reduces 

to , • 

£ 


F(2, X ) = £ 


e(Q 




-£ 


' 2 
a - 


which coincides with the corresponding expression of Section I. Cal- 
culations show that F(£, x|) changes its sign from plus to >minus 

off 

(i.e., shielding vs. amplification of when c < z Q + i for ?J h = 
1 ( see Figure 11). 
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CHAPTER III 

Higher Order Approximations for the Problem 
of Two Col inear Cracks 

1 . Linear Approximation 

In this chapter the problem of Chapter II, Section I, will be 
considered (Figure 5). All the assumptions of that section are 
assumed to hold, but instead of constant approximation of the 
resulting stress field g(x) on the microcrack, linear approximation 
will be taken. Thus, the elastic stress field on the microcrack is 
assumed to be of a form 

o(X') = o’(C’) (X' - C) + o(C’), X’ € (C - r , C’ + V) (3.1) 


where 0 '( c') is the derivative with respect to x' at the center of 
the microcrack, x' = x/£ q , c' = c / £ q , etc. represent nondimensiona) 
coordinates . 

According to the theorem on polynomial- conservation (Willis's 
theorem) [3$] a polynomial loading produces elliptic crack opening 
displacement multiplied by the polynomial of the same degree as the 
loading. Using this result the double layer potential density b(0 
can be written in the form 


b (C) = [b + b AC - C’)l 4e(Qn(U 


(3.2) 
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instead of (2.2) in the previous chapter. Unknown coefficients b Q 

and are of length units. By formula (1.3) the displacement 
vector may be represented as follows: 

C+2. 

U (x) = [b + b (£ “ C)] 4 e(0 n(€)*$(S, V d[ - 

~ ~ J ° l 

c-i 


where «(£)•*<£. x) = X) = * 2 g(S. X) 

Application of the stress operator T x to the displacement vector 
u(x) results in the stress field of the microcrack. The a^(x) com- 
ponent of it may be written as follows: 


C+'i 


o 22 (x') 


E_ 

iff. 


C-f. 


e(Q 

(V ~ x') : 


[b 


b l(C 


- c ' )]' d r. 


(3.3) 


Formula (3.3) is valid on axis (i.e. x^ = 0). In this chapter we 
will write x instead of x-| . Both coordinates x and g are on the 
horizontal axis. The integral in (3.3) diverges when x«(c - l, c + 
y.) and must be understood in the sense described in Chapter I, Sec- 
tion 2. Evaluation of this integral gives. 


°22 (X,) = °22 (C ' } (X ' " c ' } + °22 (C,) = ! [b o + 2b l (X ’ ' c ’ )] (3 * 4) 


and comparing coefficients of linear function we obtain the rela- 
tions 


°22 (C,) = X b o 


o’( c ’) - 2j bl 


(3.5) 
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Substitution of (3.5) into (3.3) expresses ^(x) as a linear 


combination of ^(c') and 22 ^ c ^ as ^°^ 0wS: 


°22 «'> ’ ° 22 (C '>t ==^ 2 

Vl '‘x-cr 


- 1] - l (C 1 - X')o’ (C')f( . - 1 ■ ■ 2 

2 22 /I _ (_L ) 2 

Tx-c' 


- i) + 


( * ( xfe )2 ■ 1)] 


(3.6) 


X £(C - 2., C + £) 

Equation (3.6) has been obtained from (3.3) by means of integration. 
(The integrals are evaluated in Appendix II.) 

Following the procedure in section 1, Chapter II boundary condi- 
tion (1.7) on the microcrack takes form (2.3) and, consequently 


« eff 

°22 (C)= - K 1 


1 


1 


/2ir(?.+G) 2(£'+6’) 


( 3 . 7 ) 


The equation for effective stress intensity factor (1.9) takes the 
form of (2.4) with o^U') given by (3.6). The system of three equa- 
tions (2.3), (2.4), and (3.7) is a system of three linear algebraic 

eff 

equations for determination of o^U), and • Substitu- 

tions of (2.3) and (3.7) into (3;6), and (3.6) into (.2.4) result in 

rff 

the linear equation for with the solution (2.14). In, this 
problem q = q Q + q^ , 


where 

f o (2,6) 

q = — ■ — , and 

° V2(r+6') 

qi . TT /2 ( H ’ +6 ’ ) 4(?/+6’) 


(3.8) 
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(both q Q and are nondimensional ) with 
1 


f(£,6) - 


| [Hi I o (c '-x , )« 


-l 


(3.9) 


f x ( 5. » 6 ) = 


[Ii(C’.X') + I 0 (C’.X’)(C’-X')]dX 


-1 


where 


C+£ 

f 


I (C'X') = £ 

O 


ds = h — i— 1) 


c-£ 


(C-X)' 


■ ( x^c r 


(3.10) 


and 


c+£ . 


I 1 (C , ,X') = £ 


e(C) 


i^dC = TT ( X ' - C’) (/l-(^) 2 - 1) (3.11) 


c-f 


The coefficients q Q and q^ , are responsible for increase of effec- 

rff 

tive stress intensity factor k^ due to the elliptic microcrack 
opening and due to the linear deviation of elliptic shape, corre- 
spondingly. 

General superposition formula (1.2) gives the resulting stress 
field in the form 


off <f>[9(x)] 

o(X) . K?" ( = — + 

/2itr (X) 


(3.12) 


4 £ 
E 


C+£ 


*^2jt(£+6) ~ X 


/ 

C-£ 


[1 - 


C-c 


4T' £+6) ] e(On(C).$(£,X)d£} 


Formula (3.12) has the, same structure as (2.15). The difference 
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appears because of the second term in brackets in the integrand. 

This term represents the correction to constant approximation due to 

the first term of the expansion of the resulting stress field on the 

microcrack (3.1). The graph of -L- for linear approximation is 

K 1 

given on Figure 13 together with the k° for constant approximation. 
"Linear approximation is represented by the upper curve. 



FIGURE 13 

The deviation of the upper curve becomes significant for small S/a, 
i.e., when the distance between the macro and microcrack becomes 
smal 1 . 

It should be noted that the solution obtained above gives low 

gf f 

estimate for both k-j and the stress field o(x) in (3.12). This is 
illustrated by Figure 14, where our solution corresponds to the 
approximation of the resulting stress field component ^(x) by the 
tangent at the center of the microcrack. 
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FIGURE H 

The results of calculations by the exact formula from [38] 

together with the upper (i.e., the approximation of the resulting 

stress field by the cord drawn through the end points of the micro- 

crack) and the lower bounds are shown in Figure 15. The formula 

€ff 0 

from [38] employed for the evaluation of the ratios in 

the coordinate system of Figure 15 may be written as follows 

**11 1— a +6+ 3Jc) (2£„ +6+ 21) - 

K° /6(6+2£) ° 

-(2 1.^+6+2 £,) 2 - (2l+£)(4 0 +t)) 
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where 


F(K) = 


dt 


A - K^sin^t 


and 


n(n,K) = 


dt 


(1 + nsin^t) A - K^sin^t 


are the complete elliptic integrals of the first and third kind, 
respectively, and 


TT 

2 


I(n,K) = 


j 

o 


dt 

(1 + nsin^t)^/! - K^sin^t 


with 


n = 


2g, 0 

X + 21 * 


and 


K 


2 V 

( 2 * 

o 



6/| 


FIGURE 15 

The lower curve in Figure 15 represents the tangent approximation, 
the dashed curve corresponds to the exact solution, and the upper 
curve represents the cord approximation. 

Thus, for the microcrack (small in comparison to the macrocrack) 
located within the range of the macrocrack asymptotic stress field 
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(2.1) the solution by means of linear approximation obtained above 
is in a good agreement with the exact solution. 

2. Quadratic and Cubic Approximations 

This section is a direct continuation of the previous one; the 
quadratic and cubic approximations of the resulting stress field 
a(x) will be considered. Our goal is to develop the method of con- 
structing the higher order approximations and possibly to represent 
the resulting stress field g(x) as an infinite series. The Tatter 
will be done in the next section. It will be shown also that under 
certain assumptions about the character of the stress field g(x) 
this series represents the exact solution to the problem. 

The elastic stress field on the microcrack is assumed to be of 
the fom« ? (x ’> s S(C’) + S*(C»)(X’-C’) + 


o" (C ' ) 


(X'-C) 
2 ! 


+ o'" (C'> 




(x’-c’y 

3! 


(3.13) 


X«(C-fc, C+2) 

where 0 '(c'), o"(c') and 0 ’ " ( c ' ) are the derivatives of the stress 
field a(x' ) at point c' . By Willis's theorem [353 the crack opening 
displacement may be represented as a third degree polynomial super- 
imposed on elliptic crack opening, so that 

2 

, b(X') = [b + b (X’-C’) + b + 

° 1 t 

(3.14) 

+ b 4 (X')n(C’) 

- 3! 

The boundary conditions (1.7) have the form (2.3) (i.e. o^U) - 

off 1 

K, ) . which leads to the conditions 

1 + 6 ) 


45 


(3.15) 


° 22 (C) 


ef f ' 

K® S (£+6), 


° 22 ( C) 


K^ f f S" (S-+6) , 


° 22 (C) 


K® ff S"' (J.+6) 


where 


s(x) = 


✓2irr(x) 


(all the derivatives s', s'', s''' are of the same units as S(x)). 
The stress field (3.13) must be expressed in terms of crack opening 
coefficients by means of (1.6). This gives for CT 22 ^ ^ ^ 

..... I (n) (X'-C’) n . 

°22 <X > * n=0 °22 <C > 


n=0 22 

C+£ 


3 K 

n(X')n(X')T x ] R | 0 b R (C - — e(C’)n(C) .*<£’ ,X’)dC 


C-x 


K! 


(3.16) 


X€(C-2-, C+2.) 

The righthand part of (3.16) has to be a third degree polynomial, 
and by comparing the coefficients on the right and on the left of 
(3.16) the relationship between 0 ^ ( 0 ) and b n may be estab- 
lished. This is linear relationship and, consequently, can be 
written as follows in variant form 

{ 0 } = (b) (b) (3.17) 

or in component form as 


(C) K=0 B nK b K 
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where {B}’ 1 is transformation matrix and { 0 } and (b) stand for 

/ \ » H •*» 

vector-columns o 22 (C): (o 22 , o 22> o 22> o 22 > 


V {b l* b 2 * b 3 * b 4 } ( 3 . 18 ) 


or, evaluating the integrals in (3.16): 




0 

2 


\ 1 ' 2 0 \ 

0 |t' 2 


.0 0 3 0 

'0 0 0 4 


(3.19) 


Note that the only two singular integrals present in (3.18) have 
been evaluated already in section 1 , see formula. (3.4). In the case 


of quadratic approximation the matrix degenerates Into 3x3 matrix 

which can be obtained from by crossing out the fourth column and the 

fourth row. It can be easily checked that for a linear approximation 

the matrix degenerates into 2x2 which can be obtained by crossing 

out the third and the fourth columns and rows (see formula (3.5)). 

eff 

In order to obtain the effective stress intensity factor k^ 
and the stress field g(x) we need to solve the system of equations 


(3.17). The inverse matrix has the form 


(B) = l 



(3.20) 


47 



The column { b} of crack opening coefficients may be represented now 


as follows 


{b} = {B} {a} 


(3.21) 


Substitution of (3.15) into (3.21) gives { b} as a linear function of 


(b) = K* ff (B) (S) 


(3.22) 


where a vector column { S > is represented by the components S(«- +6), 
S' ( £ + 6 ) , S" ( £ + 6 ) , S" ' ( £ + 6 ) . Substitution of (3.14) with { b) 
given by (3.22) into (1.6) and the latter into equation for effec- 
tive stress intensity factor (1.9) results in the linear equation 

• * v 1 t , 

ef f 

for k 1 , with th solution (2.14). In order to carry out this proce- 
dure we substitute (3.14) into (1.6) at first, and after rearrange- 
ment of some terms obtain 


22< X '> ■ - n { f e -^ ' b o + b i<*'-c'> + h. b 2«’- c ’> 2 + 


■I j. f ® (t) 


+ 57 b 3 (X , -C , ) J ]dt - ^ [b 1 + b 2 (X’-C’) +■ 


i K Cv • n imj. f _ Ai 


+ 5 b 3 (X'-C')]dt + e Ct) Ijb 2 + b 3 (X' -C 1 ) ] dt ’ 


(3.23) 


- 3 ! I t*.e(t)b 3 dt’} 
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/ f 

where t = £- x, e(t) = |1 -( — - — ) , and integrals are taken from 
(c'-t'-x*) to (c + t' - x') introducing the notation 


1 , x * ) 


ul 


1 

A ' ( x T ^c T ^ 2 


1] 


i x (r , x ’) 


n (X'-C * ) [ 1 - H 


i 2 (£\x’) 


C'+?'-X' 

e(t)dt' 

C'-i'-X’ 


tr t' 2 
2 


1 3 ( ^ * »X' ) 


C' + K-’-X' 

t ’ e ( t ) dt ’ = — j— ( x ’- c ') 

C’-JL'-X' 


(note that I and 1^ are the same integrals as in section 1, formu- 
las (3.10) and (3.11)) formula (3.23) may be rewritten as follows: 


°22 (X) --n ll o b o +t<X ’- C ' )I o +I l )b l + 

+ ttt l(X’ - c ’) 2 + 2(X’ - C)l. + l ? )b + 

25 122 (3.23a) 

+ i-[(X’ - C’) 3 I + 3<X’ - C’) 2 I. + 3(X’ - C’)I + IJ) 

J ! o 1 2 J 


Substitution of (3.22) into (3.23a) results in the following expres- 


sion 


o ' 2 < x > - - jjtys - k s "> + 5 <s ' - sW 5 "'* 


[(X’ - C')I o + Ijl + 2T3 S "l (X ' - C ^ 2l 0 + 2(X ’ " C,)I 1 + 1 2 ] 


(3.23b) 
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rH f ^ 


+ _Ls , "l(X' - C') 3 I + 3(X' - C’) 2 I. + 3(X' - C*)I 2 + (3.23b) 

3 1 A ° . 


Substituting (3.23b) into (1.9), results in 

K i ff s k ° * k f f <% + <ii + + 

pff pff off pf f 

where k-j q Q , k^ q^ , k^ q 2 , and k-j q^ represent the increments of 

rff 

effective stress intensity factor k-j , resulting from constant, 
linear, quadratic and cubic approximations, respectively. The 
solution of the last equation appears in the form (2.14) with 


q - Z. q, (3.24) 

i = o 1 

Each of q i - is defined by the structure of (1.9) and (3.23b) and 

may be written in the form 

1 


, i’2 / 1 + Y > 

•T (s - 3T4 s ) l/rrx* V r ’ x ’> dx ’ 


1 ^-o 1 £ ' 3 

a * i • — - (is* ^ c«»n 

q l IT TT K 2 ' 2*3*4 s J 


/l + X' 

x* iy*\ x,) (x ' - c '> + y*'. x 


tt 7i 2 ! 3 | ij 1 - x ' ^ X ” ^ ) y^. X') + 2 (X ' - C')^ (£', X') 


-1 


1 O S f f 1 X V I *5 *y 

V • - * m ] $ T^rr t(x ’ - c ’> l Q + 3(x ' - c ’> X 1 + 3 <x’ - c* ): 

-l 

The formulas above can be easily reduced to the ones in the linear 
approximation by setting S ' ' =S * ’ ’ =0 . 
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' ) ] dX ' 


+ I 2 ]dX’ 


2 + I 3 ]dX 



Using general superposition forumla (1.2), formula for effective 
stress intensity factor (2.14) with q^-th (i = 1,2, 3,4) given by 
expressions above, and formula (3.23b) for the resulting stress 
field o ( x) for cubic approximation may be written out. We do not 
do it here, but in the next section the general formula for the 
stress field o(x) is presented in a form of infinite series by the 
derivatives of dominating field at the microcrack center, i.e. at r 
= £ + 6 • 

In the considered case of cubic approximation of the elastic 
stress field g(x) on the microcrack all the calculations have been 
done exactly. In case of higher order approximations the complexity 
increases and only approximate calculations have been done. It 
should be noted that integrals of the type of I Q , I^ I^, etc. 
always can be evaluated in terms of elementary functions [37] and 
the procedure can be carried out in the case of- higher 
approximations also. ' ' • 

The concluding remarks of Section 1, Chapter II are applicable 
to the problem of this section also. It means that the equation for 
the effective stress intensity factor may be solved by iterative 
procedure which corresponds to the expansion of (2;14) into geo- 
metric series. Physical meaning of the terms in geometric* series 
remains the same. 

Thus, the solution to the problem under consideration can be 
refined in two directions: higher order interactions taken into 
account, and higher degree of approximation of the stress field g(x) 
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i nvol ved . 


3. The Series Solution (Two col inear Cracks) 

In this section the procedure of obtaining higher order approxi- 
mations developed in the previous section will be extended up to the 
P-th approximation, where P is an arbitrary number. It will be shown 
that for an analytic stress field o(x) the P-th approximation, 
becomes the exact solution in the limit P + <*. In the logic of 
development of this section we follow the pattern of the previous 
one. Thus, the resulting stress field g(x) on the microcrack i 
will be represented by an expression 

o 22 (x) = Z o( £ ) (C) , x«(c-£, C+£) (3.25) 

K-0 

then by Willis's theorem [35] the COD on the microcrack may be 
represented as follows 


b(X) 


= 4n(c) 


P 

Z 

K=0 


K 


(X - C) K 
K! 


(X) 


(3.26) 


Boundary conditions (1.7) and the ones obtained from them by means 
of differentiation (i.e., analogous to formula (3.15)) take the 
form: 


0 22 (C) = K® ff 


pff 

° 22 (C) " K 1 


o^CC) ■ Kf f 


S(£ + 6) 

S' (£ +6) 
S (P) (£ + 6) 


(3.27) 
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The relationship of the vector a^(t) = { 0 } to the vector 

b^ = { b } is linear and given by (3.17), where the transformation 
matrix { B } ” ^ is defined by the expression 


C+£ 

-l f 

{B} (b) = n(X)n(X)n(C) *T 

~ x J 

C O 

-X 


P K 

1 --ir r - e(D 9 (c, X)dc (3.28) 
K=0 K K ’ 


which is analogous to (3.18) of the previous section. Formula 
(3.28) (just as (3.18)) do not give the explicit form of transforma- 
tion matrix. In the previous section, all the integrals in (3.18) 
have been evaluated and the resulting structure of (B) - ^ appeared as 
(3.19). In the general case under consideration higher order inte- 
grals in ( 3 .28 ) were estimated with the help of the mean value 
theorem. The structure of transformation matrix {B}" 1 (and { B } ) re- 
mains similar to the one given by (3.19) and (3.20), respectively). 

In order to complete the scheme the equation for effective 
stress intensity factor (1.9) must be used; it is linear equation in 
view of (3.22) , which preserves its form with { B } - 1 defined by 
(3.28) and the vector { s} extended correspondingly. The solution of 
(1.9) has the usual form (2.14) with 


P 

q = I 
K=0 


q K 


(3.29) 


where 



P 

I 

i=0 


Ki 



C+£ 

r 

1 J 

c-l 


(€ - c) K 

K! 


-i. ^ -.dOdX (3.30) 
(e - X) 2 


X (C-£, C+£) 
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With effective stress intensity factor k-j determined, the last 
step left is to use general superposition formula (1.2) for determ- 
ining the stress field Q (x). 

+ 


c (x) - tcf E 

1 /2irr(X) 


P N 

+ ii 1 1 B nK S K T X 

n=0 K=0 


C+t 


C-8. 


,r r N n , (3.31) 

— e (C)n(C) *4>(C» X)dU 

n ! 


Formulas (2.14) (with q defined by (3.29) and (3.30)), and (3.31) 
are the resulting ones for P-th approximation. The rest of this 
section is dedicated to the analysis of the structure of matrix { B) 
and the conditions under which formulas (2.14), (3.30) can be ex- 
tended to exact results. 

It should be noted that evaluation of the integrals in (3.28) is 
of vital importance because only that gives the explicit structure 
of transformation, matrix . The integrals may be represented in the 
form 

C+H : ' ' ' . • . 

(X-C ) n " 2 e(Qd£, Xc(C-S., C+S-) 2 

C-K. 


I (£, X) = - fi- 
ll 


(The first four of them are Iq, 1^ , I^, of the previous section.) 

Only 1^, I-j are singular integrals and they have been evaluated 

o 

already. All the integrals I are of the type R(x, x + ax + e) 
where R is rational function and, therefore [37], can be integrated 
in terms of elementary functions. It is enough for our purposes, 
however, to estimate them by the mean value theorem 


e (E ) ' , , n-i 

!„(*» x) = n - 1 {[X " (C + £)] - [X - (C - &)] } 
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where e(x Q ) _< 1, for X Q €(C - 2 , c + 2 ). The last formula can be 
rewritten as follows 


2 (X ) n-1 

i ( 2 , x) = - r ~ r 1 

n n 1 tn=0 




c- ,<X-C) n ' m - I n ,n [(-1)"-Il 
n-1 


(3.33) 


where C m , are the binominal coefficients, 
n- 1 

Using the expression (3.33) with the exact results for the 
integrals I Q - I 3 (section 2) the stress field on the microcrack can 
be expressed as follows: 


N 


°22 (X) = 


- n 1 E 0 <-> It c : (x - c) " + , 

n=0 n=l 


N 

Z 


nC 1 ^ (X-C) n + 
n n ! 


(3.34) 


+(-l) 


K+l 


N 

Z 

n=K 


C K bn 
n n ! 


e (X ) 
o 

K - 1 


K-l 

+ Z 

m=0 


C K-1 ^ ( (- 1 ) m - 1 ] (X-C) n 


By comparison of the coefficients at the same powers of (x - c) in 
the righthand part of (3.34) and its lefthand part in the form 
(3.25), we obtain the transformation matrix (B}~\ from which matrix 
(B) can be obtained. 

The structure of both matrixes {B}” 1 and { B } appear to be the 
same as in (3.19), (3.20), i.e.. 


-1 


{ 8 } 


. E 
'• 2 


/I 0 “o? 1 0 o 0<. 1 * 0 °06 l b 

0 2 0 Ojjl 7 0 Ojjl 4 0 

oo 3 o 7 o o, 6 i 


A 


(3.35) 


\0 0 0 


./ 
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and 


/ 


{B} = £ 


0 •oi* 


\ 0 0 


e o/* 0 V 


B n* 


•2 


•t 


't 


0 B »* 


h/' 0 *V’ 


\ 


:/ 


(3.36) 


where coefficients ^ , » Q4 , <> 06 , etc. and 8^, 8 Q4 , 6^, 
^13’ ^15’ etC * Can ca ^ cu ^ atec * ^ rom (3.34). It should be noted 
that coefficients ij decrease by both indexes i and j. 

The structure of transformation matrixes (B) - ^ and (B) described 
above is not a result of approximate calcultions of integrals I n « 
This can be shown by considering the integrals I n : Notice that the 
product of binomial expansion of the term (x - £ ) ' multiplied by 
elliptic COD e(^) represents the sequence of odd and even terms. 
Integration within symmetric limits leaves only even terms. 

Formulas (3.35) and (3.36) complete the consideration of the 
matrixes (B}~^ and { B }. If the o^Ox) component of a stress field 
g(x) is analytic in the interval (c - i , c + t), then (3.25) may be 
considered as its Taylor's expansion when P -*•<» . Correspondingly 
(3.26) is Taylor’s expansion for b(x) and formulas (2.14) (with 


given by (3.30)) and (3.31) represent the exact solution of the 
problem if the series in (3.30) and (3.31) converge. 

The convergence in both formulas is defined by the behavior of 
the following series 
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f(£.x)dC 


(3.37) 


Z Z 
r>>k k>0 


c+£ 

' 

B . S. 
nk k 

c-£ 


a-c) n 

n 


where f(c , *) = -r> is a scalar function for formula (3.29) and 

U-xr 

f(S, x) = e(c ) n(C).4>U, x) , where x*(c - £, c + 2 ,) is a vector 
function for formula (3.31). Taking into account that 


<-1 > k (2k-d i : 

— — . - f 

r = £+6 /2 tt ( £+6). ( £+6)^2^ 
for K>0 

and that 



(3.38) 


c+£ c+£ 

r . n . n 1 r 



r . , x n n 

U -•7- c) ~ f«. x)d£ | < ^ 


► 

f ( C . x)d£ 


n. - - 1 — n. 

t 



(3.39) 


(3.37) can be bounded as follows: 


c+£ 


Z B . S. 

n> k-0 nk k 


— -~t — fU. x) 


c-£ 


c+£ 


(3.40) 


Z Z B 


(2k - d:: 


n b b 

✓2tt(£+ 6) n>k>0 2 n! 


(1 + 6 /£) 


K 


*U. x)dC 


c-£ 


The series on the right of (3.40) is absolutely convergent and the 
rate of convergence depends on 6/z : the bigger the ratio <5 /£, the 
better the ' convergence is. This result is understandable from a 
physical point of view; the small ratio indicates that either the 
microcrack is close to the main crack tip (i.e., 6 ->• 0 ) , or the 
microcrack length £ is large; . in both cases many terms in the 
stress field expansion are needed. For a large ratio either the. 
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microcrack is far from the main crack tip (i.e., s is large), or the 
microcrack is small ( i .e. , 0) ; in both cases the correction to the 

main crack dominating field is small and our series (3.37) converges 
rapidly. 

As conclusion to this section we will formulate the above as the 
following statement: if a system of two cracks in linear elastic 

medium under the assumptions of plane problem and mode I loading 
located as on Figure 5, and the microcracks are imbedded in the 
dominant field of the main crack, then the analytic stress field 
g(x) can be represented as an infinite series (3.31) where: 1. the 
effective stress intensity factor is given by (2.14), (3.29) 
and (3.30), 2. the first term represents the dominating asymptotic 
stress field of the main crack, 3. the second term represents the 
series by the derivatives of the dominant field taken at the 
center of the microcrack with coefficients which depend on small 
crack length and location, 4. the series is absolutely convergent 
and the rate of convergence depends on the ratio of the distance 
between two adjacent crack tips 6 to the smaller crack length z . 

The obtained result is an exact series solution to the problem 
of this chapter. It should be noted, that the exact solution to 
this problem has been obtained in [38] by means of complex variable 
techniques. It will be shown in the next chapter that the method 
suggested in this section may be generalized and the exact solution 
to the problem of a macrocrack surrounded by an array of arbitrarily 
located microcracks of arbitrary lengths and orientations can be 
obtained. 
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CHAPTER IV 


The Multiple Crack Interaction Problem 
1 . Stress Field in a Neighborhood of a Macrocrack Tip Surrounded by 
an Array of Microcracks 

In this section we consider a plane elastostatic problem of the 
interaction of a macrocrack under mode I tensile loading with an 
array of rectilinear microcracks at the macrocrack tip. 

Let us introduce the following notation. 


1 . 


o (x) = K 

O l 


$[eqp] 

/2iTr (X) 


(4.1) 


is the asymptotic stress field of a macrocrack of length 2 £ q with K-j 

- as stress intensity factor 4>[9(x)] - as asymptotic angle distibu- 
tion tensor, and r(x) - as a distance from macrocrack tip to the 

. t.W*)* .r>0t> tln) <x*) 


point x. 

2 . 


(4.2) 


is m-th derivative of traction vector t in the direction of s-th 
microcrack evaluated at the center of the microcrack. The direction- 
al derivative of higher order must be undestood in the following 
sense: if - X^q * a(x* - x*g) is the equation of the rectilinear 
segment of s-th microcrack (see Figure 16), then the stress g(x*, 
xp on that segment can be represented as o[X^, X 20 + a(x l X 1(P^ 

- i.e. a function of x* coordinate only. Then the* directional 
derivative takes the form 


a (n,) ( x s ) - -L 
~° dx'J 


°[x*. 


20 


“( X i - x iV ] 


(4.3) 
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FIGURE 16 

3. b£ n) - (b) 

stands for the set of vector coefficients of a polynomial which is 

k 

superimposed on the elliptic crack opening of fc-th microcrack e( £; ) 
(see formula (4.14)) . 

4 - <?> - B «, (X 0- X o> < 4 - 4 > 

is a linear operator which depends on the positions of centers of 
s k 

two microcracks x , x . and characterizes the influence m-th deriva- 
o o 

tive of s-th microcrack stress field on n-th derivative of k-th 
microcrack opening displacement coefficient (b). The linear 
operator {B}is a second rank tensor (subscripts in (4.2) and in all 
the operators below are not tensorial). The definition of the 
linear operator { B) is given below. 
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(A(X)} = A (X K ) 
~ ~ ~r> ~ 


T 

~X 


(C K - 


/r K x W -K 
e(C ) 4> ( t, , 


X)d£ 


K (4.5) 


K 


n ! 


is a matrix which characterizes the stress associated with n-th 
derivative of the microcrack opening coefficient {b} of k-th micro- 
crack at arbitrary point x. Matrix fA(x)} is the matrix of the micro- 
crack arrays In formula (4.5) T x - is stress operator transforming 
the displacement field u(x) into stress field o(x), e(^) - is a unit 
elliptic crack opening displacement, i>(§, x) - is the second Green's 
tensor defined in Appendix I , £ is a coordinate on k-th microcrack 
of length 2fc k ; matrix {A( x)} is a third rank tensor, 

6. U°} = A°(X K ) (4.6) 

is a linear operator which characterizes the increase of stress in- 
tensity factor K° due to stresses associated with microcrack opening 
coefficients {b } of the array. The linear operator {A 0 } is a vector. 
The definition of {A 0 } is given below. 

Using the notation above, the following statement can be proved: 


If, 


1. The microcrack array has been located in a close vicinity of 
a macrocrack tip and the characteristic linear dimension of the 
array is small in comparison to the macrocrrack length (i.e., small 
scale microcracking model ) , and ; - 

2. The resulting stress field o(x) induced by the macrocrack 
and the rest of microcracks on the line of each microcrack can be 
approximated by the polynomial function. 

Then the resulting stress field o(x) can be fully characterized 
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by the asymptotic stress field of the macrocrack ? 0 (x)» the values 
of its derivatives oj m ^(x o ) in the directions of microcracks 
evaluated at the centers of microcracks, and the second Green's 
tensor <j>(g, x) as follows: 

o(X) = ({1} + {B} { t } {A°}.)o o (X) + {B} { t > {A(X) } ( 4 . 7 ) 


or in index notation 


?QP = [1 + E B m (x* xK H (x s )a°(x K )]o (X) + 

~ ~ m,n,s,K mn ~o ~o 'in ~ ~ ~ ~o ~ 


■ (4.7a) 


m,n^s,K 5 m «o- So>£„«o> V* K > 


Proof 

The elastic stress field of the system of cracks under consi- 
deration can be represented by means of general superposition 
formula 


a(X) = a(X) + . 1 . o. (X) 

~ ~ - ~ 1=1 ~i ~ 


(4.8) 


where 


, , ef f ^ 0( ^ )1 

o(X) = K 

1 /2Tir(X) 


(4.9) 


and 
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(4.10) 


o.(X) = T x | b.(Q- *a, X)dC- 

a. 

i 

with as a double layer potential density which can be inter- 

preted as a crack opening displacement. Fprmula (4.9) represents 
dominant stress field of the main crack with - as effective 
stress intensity factor which reflects the impact of the microcrack 
array on the macrocrack. It appears as an unknown quantity in (4.9) 
along with unknown components of vectors b^(§), where i = 1 , 2 , 
... N. 

Equilibrium equations are automatically satisfied for the re- 
sulting stress field a(x) in the form (4.8). The stress field of 
N 

microcracks g.(x), which is defined by (4.10), satisfies the 

equilibrium equations because of the properties of the second 
Green's tensor $( 4 , x) . The dominant stress field 6 (x) , given by 
(4.9) satisfies the equil ibrium equations because of the properties 
of the asymptotic crack tip solution. 

The equations to be satisfied are the boundary conditions on the 
macrocrack and each of the microcracks. For the small scale model, 
the boundary condition on a macrocrack can be substituted by the 

rff 

equation foh effective stress intensity factor K-j which fully 
determines the asymptotic stress field of the macrocrack. 


0 fe + x 


K eff = + JL 

1 M 


— nWnW.^g.QOdx (4.11) 


0 -P- 


Boundary conditions on microcracks appear to form a system of 2N 
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singular integral equations for 2N unknown components of N vectors 
of double layer potential densities b(x) (K = 1, 2, ... N) . 


?i { ° ( V + K =i I x k K <§>*!(§. ?i )dC + 

o 

X, 

b (§)• 1' (§, X )dC) = ° 


(4.12) 


where x. - is a coordinate on i-th microcrack. Equations (4.11) , 

pf f 

(4.12) constitute a system for determination of 2N +1 unknowns 
and the components of N vectors b^E ) • 

In order to solve the system of equations (4.11), (4.12) let us 

represent the polynomial stress field o(x) on k-th microcrack by: 


P rx K 

o(x K ) = I 0 (n) (x K ) - £_ 

~ - n=0 ~o o n! 


(4.13) 


Formula (4.13) and Willis’s theorem [35] permit to represent the 
1/ 

GOO b(x ) as a following polynomial 

P 


(X K - x K ) 

. _ I An) V 

5 (X 5 - n=0 ~K ( ?o } nl 6 (X) 


(4.14) 


In view of (4.13) and (4.14) the boundary conditions (4.12) may be 
written as follows: 


n(X K )T v 


P N 


P (T K - f K ) n 

i t/ n ' a„) — e(c K ) * a,, xi = 

n=0 ~° ~ ' K 




,K, 


-I l n(X ') • T 

, ~ -O ~X 
n=0 s=l 


S (£ S - O s s K (4.15) 

b <n) ( So } „ i ~~ e(C S ) * (C S , X K )df. - 


- K® ff S(X K ) 


where S(X K ) = n(X K ) ^11 

/27ir(X) 
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Expression on the left of (4.15) represents traction on K-th micro- 
crack which, by assumption, may be represented as a polynomial . 
(The lefthand part of (4.15) contains singular integral, when n=0, 
which converges in Cauchy's sense). Consequently, the righthand 
part of (4.15) (i.e., the traction on K-th microcrack induced by the 
rest of the microcracks and the dominant field) is a polynomial 
also; then (4.15) can be rewritten in the form of a system of 
equations for the coefficients of the polynomials of the lefthand 
and righthand parts of (4.15). The procedure described above can be 
carried out by differentiation of (4.15) at point x q in the 
direction of k-th microcrack P times. The differentiation: results 


in the following equation: 


„£o A !«<?:. £$ * - & M 9„X- *„>!{ 


K l ff VA 


(4.16) 


where S m (x^ is m-th derivative on s-th microcrack, m = 0 , 1 , 2 , 
... P, S = 1, 2, ... N. b" = (x^) . The 1 inear operator H n 


^ c U 

(XqXq) is given by the following expressions: 


H (X s , A =| T . 

Mnn ” O O ] . v 


when s = k 


a s - c s ) n 


e(C s ) 4> (m) (C S , X K ).n(X K )d£ S (W'i)' 


taken at point x = x. 


when s # k 


where ( 5 s , x^) is m-th directional derivative of Green's tensor 
at point k taken in direction of k-th microcrack defined by (4.4). 
Similarly, the linear operator I mn (* 0 » * Q ) is given by the expres- 
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Sion : 


I (x s 
-mn 'Co ’ V 


r. 


T- 

~X 


' i 


tt s - C B ) n 


n ! 


when s / k 


e(§ S )> (m) (§ S , /).n(/)d^ S (4.18) 


when s = k 

The set of elements in (4.16) constitutes a P x N matrix. In 
symbolic notation equation (4.16) may be written as follows: 


{1} - (b) = - { H } *{b} - K^ ff - (s) 


(4.16a) 


eff 

together with equation (4.11) for stress intensity factor k-j boun- 
dary conditions on microcracks (4.12) form a system of (P x N + 1) 

linear algebraic equations with P x N unknown components of matrix 

eff 

{ b} and unknown quantity . Substitution of formula (4.10) with 
b^(c) given by (4.14) into (4.11) results in . 

K t" ■ K l + (A °> • (b > (4.19) 


where linear operator {A 0 } is given by the expression 


{A°} 


Z 

r° 


/u£ 


Z + X 
o 


to -l 


Z - X 
o 


A 

n(X) s E . H (X S . X)dx (4.20) 
' - n=l s=o -on - ~ 


with H on (x^ defined by (4.17). Formula (4.20) defines linear, 
operator {A 0 } of (4.5). Operator (A 0 ) characterizes the impact of a 
microcrack array on the main crack. 

Substitution of (4.19) into (4.16a) gives a system of linear 
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algebraic equations for the determination of (b) 


{1} {b} = - ( {H} + {§} 4°} ) {b> - { t) 


(4.21) 


The last system of equations yields the obvious solution 

(b) = {B} { t } 


(4.22) 


where 


{B}=-({I} + {H}+{S} {A°}) -1 (4.23) 

Formula (4.23) defines the linear operator {B } of (4.4). 

In index notation (4.20) takes the form 



N 

I 

s=l 


? s K s 

1 B (X , X ) t (X ) 
m=o -vmn-o -o ~m 'O 


(4.22a) 


Formula (4.22) represents { b} as a linear function of directional 
derivative^ of the asymptotic stress field o Q (x) (n. g Q (x) = K° 

s(x). Substitution of (4.22) into (4.14), (4.14) into (4.10), and 
(4.19) into (4.9), with subsequent substitution of (4.9) and (4.10) 
into (4.8) furnishes (4.7). This completes the proof. 

Thus, taking into account the analytic character of the sol ution 
of a plane problem of el astostatics [33], the resulting stress field 
o(x) can be approximated by a polynomial as closely as desired. Con- 
sequently the solution (4.7) obtained above can be made as close as 
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desired to the exact one. 

It should be noted that each term in formula (4.7), for the 
resulting stress field 0 (x) can be given a clear physical interpre- 
tation; the first term represents dominating stress field of the 
main crack c(x) and the second term represents the stress field of 
the microcrack array embedded into the stress field of the macro- 
crack a (x). Dominating stress field 5(x), in turn, consists of two 
'o' ~ - 

terms, first of them being the asymptotic stress field of the main 
crack o q (x), and the second term results from the impact of the 
microcrack array (embedded into asymptotic stress field) on the 
macrocrack. 

The other remark concerns the method of solution of a system of 
equations (4.21). Actual construction of operator (B) may present 
considerable difficulties for an extensive microcrack array. This 
is one of the reasons for suggesting an iterative procedure as an 
alternative to exact solution. Another reason is that an iterative 
procedure has clear physical meaning and, thus, gives an insight 
into the nature of solution (4.7), as will be shown below. 

Iterative process for equation (4.21) can be arranged by multi- 
plying it by* (If* operator. This gives 

{b} (n+1) = - {I }' 1 [(H) + { S ) { A° } ] (b} (n) - {I }" 1 (t) (4 ' 23) 

where superscripts (n + 1) and (n) refer to the corresponding steps 
of iteration. Choosing {b| 0 ^ = -{1}^ {t} (i.e. microcrack opening 

coefficient { bf which results from the main crack field only) as 
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zero approximation, we obtain 

{b) (n) = [U) - {D} + {D} 2 - {D} 3 + ... + (-l) n {D} n ]{l} -1 (t) (4.24) 

where {1} is a unit operator, and 

{D} = {I} -1 [{H} + {S}{A°}] (4>25) 

The multiplication of linear operators must be understood in the 
following sense: e.g.. 


{I}" 1 {H} = 2 I -1 (X s , x K ) h (X s x q ) 
m,s -ms 'O 'O ~mp -o ~o 


(4.26) 


Consequently, double iterated operator {D? is represented by formula 


D 2 (X* 
-np ~o 


* 


I 

m,s 


D 

-mn 



X K ) D 

'O ~mp 



(4.27) 


Substitution of (4.26) (with {b)° = - {I} -1 (t) into (4.7) after 
rearrangements can be written in the form 

o(X) = O o (X) -(b) (o) (A(X) ) + {pHb) (o) {A(X)} - {p} 2 {b}(°){ A(X) } + 


(4.28) 


.. - (b }° (A}° o_ (X) + {D}{b} (0) {A}° a o (X) - {D} 2 (b) (o) {A°}a„ (X) + ... 


~o 


~o ~ 


Each consequent term in the last formula can be interpreted as 
follows: zero order term - g Q (x) , i.e., asymptotic stress field of 
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a macrocrack; first order term - {b}^ { A(x)} , i.e. stress field of 
the microcrack array (characterized by matrix { A( x) } ) embedded into 
asymptotic field ? Q (x) (note that (b)^ represents the set of 
microcrack opening coefficients for the array embedded into asymp- 
totic field o o (x)); second order terms are characterized by triple 
products: { DHb}^ { A(x)} and {b}^ { A 0 }o 0 (x); the first term 
with the once^iterated operator {D} gives the stress field generated 
by the microcrack array with microcracks subjected to stresses 
induced by other microcracks and the main crack, the microcrack and 
the main crack being embedded into asymptotic field; the second term 
with operator { A 0 } gives the stress field of the main crack subject- 
ed to the stresses induced by the microcrack array, the latter 
embedded into asymptotic stress field of the main crack. 

Thus, double products in (4.28) account for I-st order interac- 
tions, triple product - for 2nd order interaction, etc. The coordi- 
nates above can be illustrated by the following diagram. 



FIGURE 17 
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Thus, the stress field solution o(x), given by (4.7), permits 
two different types of approximations: 

1. Approximation by number of derivatives of the asymptotic stress 
field taken into account and 2. Approximation by number of physical 
interactions between microcracks taken into account. 

2. Evalution of the J-Integral 

In this section we evaluate the energy release rate J per unit 
crack layer extension, for the CL it can be expresed in the form of 
path- independent integral around the CL active zone [25]. 


J K I (f6 Kj " °ij U i,K )n j dr 


(4.29) 


The active and inert zones and Vp the line of their separation 
r(t) (i.e. the crack layer trailing edge [25]), and the contour r 
are presented in Figure 18. 


r 


FIGURE 18 




71 



It is shown in [25] that J-integral given by'(4.29) does not depend 
on the contour of integration r if the latter connects the end 
points of the trailing edge leaving the active zone inside. 
The path- independent property of energy release rate J follows from 
the principle of minimum of strain energy with the additional condi- 
tion of invariance of the strain energy with respect to translation 
of a CL. Consequently, the J-integral is path independent for homo- 
geneous medium only. For a model of the CL described above (i.e., a 
macrocrack surrounded by an array of microcracks) Figure 18 must be 
substituted by Figure 19. 



FIGURE 19 

The energy release rate J preserves its path independent 
property for the model if evaluated on the contour surrounding the 
active zone. 

Let us consider the problem of interaction of two cracks located 
on one line (i.e.. Figure 5) in order to indicate in a simple 
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problem the particular path which leads to the solution of a general 
problem. The convenient path consists of the union of the three 
contours, r Q , , and ^ and the rectilinear segments connecting 
them, depicted in Figure 20. For the loops r , , r «, in Figure 20 



the integrals on rectilinear segments cancel and only the integrals 

on r , r 1 , r 2 , must be evaluated. O-integral evaluated on the loop 

around the crack tip is shown to be [30]: 

K 2 

T _ 1 (4.30) 

E 

where Kj is stress intensity factor, thus, each of the integrals on 
r o , , r 2 , can be calculated by means of (4.30). 

This line of reasoning can be applied to the general problem of 
interaction of a macrocrack with a microcrack array (Figure 19). 
The energy release rate J for the general problem can be represented 
as follows: 

j = j + AJ ‘ (4.31) 

O 

where 
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J = 


( K ® ff ) 2 


(4.32) 


is the energy release rate associated with the main crack and evalu- 


ated on the loop around the tip of the main crack r , (it should be 
noted that the microcrack array al socontributes to d Q through the 


pf f 

use of ) and 


N 

at- E A J . 

AJ ' i=l 


(4.33) 


is the energy release rate associated with the microcrack array. 
The parts of each Aj^ associated with two ends of the microcrack 
(i.e. on the countous similar to and r^) enter with opposite 
signes (e.g., on the countour with (-) sign, and on the countour 
r~ with (+) sign) because of -the directions of the normals n. to the 

^ J 

contours. Fromt his remark follows that for piecewise constant 
appoximation of the resulting stress field o(x), each vanishes 
because the stress intensity factors on both ends of each microcrack 
are equal . 

In the general case of the resulting stress field a(x) given by 
(4.7) the energy release rate associated with i-th microcrack of the 
array can be represented as follows: 


1 2 
A J. = i [(K 
E v 10 


- 4) 


+ (K 


20 


- K n)] 


(4.34) 


where and K^q - are the mode I and mode II stress intensity 
factors at the microcrack tip which is closer to the main crack, 
respectively, and K^-j - are the mode I and mode II stress inten- 
sity factors at the microcrack tip which is farther from the main 
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crack tip. The second term in (4.34) appears because of the 
presence of the mode II loading on the microcracks [36]. 

In order to give the nontrivial example formula (4.34) has been 
used for evaluation of J-integral in the case of linear approxima- 
tion of the resulting stress field o(x) on each microcrack. 

The stress intensity factor for a crack of length a loaded 
by the normal stress o^x) defined by the following expression: 


K l = 


a 

ir 


1 + i “22 <X >‘ 1X 


1 - 


(4.35) 


-1 


Analogous expression is valid for the loading by tangential stress 

°2i ( *) • 

Under the assumption of linear approximation of the resulting 
stress field o(x) both of the above mentioned components are linear 
functions, and (4.35) can be rewritten as follows: 




(aX + B)dX 


(4.36) 


where a and b are the constants of stress units, and x is the non- 
dimensional coordinate. 

Integration yields 

K 1 * + 8) (4.37) 


The same formula is valid for but a and g would represent the 
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coefficients of linear function o^(x). 

Using the notation of (4.34), formula (4.37) gives K^, while 
K-jg is given by 

1 

£ f 

K 10 = T? 

* 

-l 

Thus, aJ can be represented as follows: 

AJ = (^a + 6) 2 - (-ia + 6) 2 ] = (4.39) 

Using the results of the previous section, the coefficients 
a and B can be determined in terms of directional derivatives of the 
resulting stress field a(x). formula (4.21), where ( b } given by 
(4.22) substituted into the righthand part, determines the traction 
and its derivatives at the centers of each microcrack of the array, 
i.e., 

(d-o}= [(H) + {SKA 0 }] { B } { t } - (t) (4.40) 


/ 


(ctX + B)dX = (- |cx + B) (4.38) 


where n is unit normal vector at the center of each microcrack. In 
(4.40) the matrix {n*a} is determined by the values of dominating 
field traction and its directional derivatives at the center of 
microcracks {t} . 

The coefficients = o'^x 1 *) and ) can t> e obtained 
from (4.40) (index k refers to the k-th microcrack). The same is 
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true for the shear mode. Substituting a k and into (4.39) we 
obtain AJ in a form: 



(4.41) 


where the derivatives of the stress components are taken with re- 
spect to dimensional coordinate, and as a result in (4.41) appears 
the multiplier £ 2 . Substitution of (4.41) into (4.31) gives the 
final formula for the energy release rate in the general problem. 

In the particular problem of macro and microcrack located on one 
line (Figure 5) using (4.41) the expression for energy release rate 
can be written as follows: 



(K * ff ) 2 



[1 



(4.42) 


The expression above is approximated for 6/)t>0.25. For smaller 6 h 
higher order approximations have to be taken into account because of 
the strongly gradient field in the vicinity of the macrocrack tip. 
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CONCLUSION 


1. The elastic stress field solution to the problem of interaction 
of a macrocrack with an array of arbitrarily located and oriented 
microcracks (within the limits of a smallscale model) has been 
obtained. 

2. The elastic stress field solution for two and three crack inter- 
action problems of interest to fracture mechanics has been obtained. 

3. It has been shown that an array of microcracks in a main crack 
tip vicinity can either amplify the effective stress intensity 
factor or reduce it depending on the array's configuration. 

4. The method leading to the solution of the general problem de- 
scribed above (i.e., the macrocrack-microcrack array interaction 
problem) refines the one suggested in the work [32]. 

5. Using the obtained elastic solution the energy release rate 
associated with the crack-layer translational motion has been evalu- 
ated . 
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APPENDIX I 

Derivation of the Second Green's Tensor for 
Plane Problem of El astostatics 

In this derivation of 2D Green's tensor we follow the routine of 
work [33] for obtaining of the second Green's tensor in 3D problem 
of elasticity. 

The second Green's tensor can be defined by the following ex- 
pression: 

5 X • g -■♦(x, §) • Q (1) 

where Q - is a unit force vector applied at point § of an infinite 
elastic plane, o - stress tensor in the plane, and n x - is a unit 
normal vector at point x of the plane. Thus, influence tensor"' <t< x ,§) 
has been defined as a linear operator which transforms the unit 

force Q applied at some point into a traction vector at the other 

- ? 

point". Surrounding the point with an imaginary closed contour we 
can write the equilibrium equation for the interior region: 

r 

| n • g dT + Q = O (2) 

r 

The value of integral in (2) does not depend on the contour of 
integration r , which may be chosen to be a circle. ' Then equation 
(2) taks the form. 
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I 5 * g dr + 9 = 0 (3) 

r 

where r is a circle of unit radius where R = x - s, \ R | = R. It 
follows from (3) that the magnitude of resultant force on any 
contour with a point inside of it does not depend on R and always 
equals - Q.. This is possible only if components of stress tensor 0 
decrease as R~^ . Consequently, the components of the displacement 
vector u must have logarithmic singularity at point §. Therefore, 
in the Papkovich-Neuber representation for displacement vector u, 

u ■ 4(1 - v) B - grad(R • B) , (4) 

vector B has to be chosen as 

B = a £n R Q (5) 

where a is an unknown scalar constant, which can be determined from 
(2) after construction of a stress tensor o on the basis of (4). 
The introduction of the harmonic scalar into the Papkovich-Neuber 
representation is unnecessary in this problem. After substitution 
of (5) into (4) the displacement vector u takes the form, 

Q • R 

u = a[{3 - 4v) Qtn R - ~ R (6) 

R 2 

The strain tensor e = 1/2 sym vu , (where "sym" symbol refers to the 
operator of symmetrization, and v operator is referred to 2D space) 
may be represented as follows: 

QR + RQ Q • R Q • R 

£ .= a[ (1 - 2v) C — ~ E + 2- — r~RR] ( 7 ) 

R 2 R R 
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where E is the unit second rank tensor in 2D space. Using the 
strain tensor (7) and Hook's law for plane problem one can obtain 
stress tensor in a form: 


- [U - 


2v) (QR + RQ - Q • RE) + 2- 


.2 •- 


RR] 


( 8 ) 


Then, equation (2) can be rewritten as follows: 


2yct f 



[ (1 - 2v) (n • QR + n • RQ - nQ • R) + 2- 


R 


-n*RR]dr* + Q = 0 
(9) 


★ 

Taking into account that R n = R on the circle r the following re- 
lations hold 


Q . R 

n . QR - nQ . R = 0, and n.RR = Q.Rn (ig) 

In view of (10), (9) takes form. 


2ua 

2 


[(1 - 2v)RQ + 2Q • Rn]dT* + Q = Q (g a ) 


r* 


Both integrals in (9a) yield 


RQ dr* = 


0 

Q • Rn dT*= 2 ttR Q 


r* 


Formula (9a) with integrals given by (10) results in 


( 10 ) 


a = - 


8utt( 1 - v) 


01 ) 
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Substituting (11) into (6) we obtain 


1 R 

U = ± [(3 - 4v)tnR • Q - — R • Q] 

8 ptt Cl - v) R 

From the last equation, it follows that the first Green's tensor 

✓ 

(i.e., Kel vin-Somigl iana tensor) defined by the relation u = U.Q can 

be represented as follows: 

1 RR (12) 

U = - [(3 - 4v) £nR • E - ^] K 1 

8 7T vj ( 1 - v) R 


Substitution of (11) into (8) gives traction vector at point x in a 
form. 




V? = 4tt ' (i vW [(1 ' 2v)( ?x?‘9 - * Q1 

R 


(13) 


From the last formula the second Green's tensor can be obtained as 


follows: * 

?<§■ ?> ■ «, R 2 ( i - u) ' (1 - 2v)< 0x? - 5?x 


n x RE) - 2 



• R 


(14) 


It should be noted that 4 >(x,s) defined by (14) can be used in the 
representation of displacement vector u in terms of double layer 
potential as follows: 

« 

u(X) = [ b(§) • X)dC 


where and x change places and b(x) is double layer potential 
density [33]. 
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APPENDIX II 

Evaluation of The Integrals 


1 . 


c+2 


I = 
o 


/TFT 


c-i 


C + 2c£ - (c 2 - 2 2 ) 
U -X) 2 


dC, X£(c - 2, c + 2) 


Substituting t = £- x, the integral takes the form 
X-(c+2) 


I = 
o 


A 2 - (x-t-c ) 2 


/r 

dt = - {y + 


X-(c+2) 

r 


X-(C+2) 


X-(c-2) 


X-(c-2) 


— + (X-c) 

A 


d t , 
t/T 


x-(c- 2 ) 


where T = £ - (x - c - t)^ 


The integral I in the last expression may be found in [37]; the 
result of integration is 


I = 11 [ X-c 

° * Ax - c ) 2 - A 


i] 


2. Analogously 
C+2 


*1 = 


C-i. 


/- p + 2c£ - (c 2 - 2^) 

U - X) 


x-(c+2) 

r A 2 - (X-t-c) 2 


dC = - 

x-(c-2) 


dt = tr (X-c) [ 1 


- 1 


X~(c-2 , c+2) 


3. The integral represented by formula (2.30) Chapter II, can be 
estimated using the mean value theorem as follows - 
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3xp - « - x*r - 6(s - 


e(£’) 


xp 2 X ’ 2 


-£ 


• 2 '23 
(<£ - Xp + X 2 ) 


dC 


e(C’) 


o • 4 

f 3X 2 


J 

-£ 


* 4 »o*o 

Xj) - 6« - Xp^X^ 


* o ' n 5“ 

IU - x 1 ) / + x 2 ] J 


dC , 


where ^ € (-a. 0 


Each of three integrals 

£ 

dt 

2 '2 3 

_{ t(c - xp z + x 2 r 

£ (£ - x') 2 d4 

f 1 

i 5 1 ? ^ 

J [(C - xp^ + x 2 ] J 


£ 

' 

< 

- £ 


U - x ’) 4 d^ 

- r~2 ' 2 3 

[U - xp + x 2 ] 



are tabulated in [37] and the integration results in formula (2.34) 
Chapter II. . 


4. Substitution of formula (2.36) into (2.33), Chapter II, leads to 
the integral 


IU + x ' , ^ 

. V” - x *£ -X- £ " £ -X+£ ^ dX 

— £ " O O O 
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where &' = n/n Q , by means of nondimensional i zation and subsequent 

2 

substitute 1 - X = 2 sin t. The last two integrals are tabulated in 
[37]. 
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APPENDIX III 


Evaluation of The Influence Function FU, x) 


The influence function FU, x) is given by formula (2.29): 

r 


F(£, X) = Vn(X)n(X)T 
- L - - - - -x 


eU)n(0*oU)‘*(5, X)d£ 


The operator n(x) n(x).J x applied to the displacement field produces 
component of the stress field. In order to evaluate the o ^ 


component of the stress field by means of (2.6), formulas (2.7) and 
(2.8) must be used. Taking into account that 


nU)-oUWU, X) = 




!iK (5 ' 


X) = 0 


$ 

22 2K 


(because n^ U ) = 0, r^U) = K and <^2 = 0 f° r this P roblem ) and 
substituting the last expression into (2.7), with <i>^ given by 
(2.8), the FU, x) with the help of (2.6) takes the form 


+ 


v 

l-2v 


F(£, X) = Z 


eCO { - — u [v(X 
R 


• 4 x 2 [x 2 + (x i - ° 2] - ^“r 

R R 


o 2 + 


v(3 + ,2v) 2 

1 - 2v A 2 


4(1+ v) 



X 2 + 



where „ ? ? 

K = (X 2 - O + 

The last equation reduces to (2.3) by factoring out the 1/R^. 
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